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Spin flips I: Evolution of the angular momentum orientation of 
Milky Way-mass dark matter haloes 
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ABSTRACT 

During the growth of a cold dark matter halo, the direction of its spin can undergo rapid 
changes. These could disrupt or even destroy a stellar disc forming in the halo, possibly 
resulting in the generation of a bulge or spheroid. We investigate the frequency of signifi- 
cant changes in the orientation of the angular momentum vector of dark matter haloes ("spin 
flips"), and their degree of correlation with mergers. We focus on haloes of mass similar to 
that of the Milky Way (MW) halo at redshift z = (logio M/ /i^'M© = 12.0 -> 12.5) and con- 
sider flips in the spin of the whole halo or just its inner parts. We find that a greater fraction of 
major mergers are associated with large spin flips than minor mergers. However, since major 
mergers are rare, the vast majority (93%) of large whole-halo spin flips (6 > 45°) coincide 
with small mass changes, not major mergers. The spin vector of the inner halo experiences 
much more frequent flips than the halo as a whole. Over their entire lifetimes (i.e. after a halo 
acquires half of its final mass), over 10% of halos experience a flip of at least 45° in the spin 
of the entire halo and nearly 60 percent experience a flip this large in the inner halo. These 
numbers are reduced to 9 percent for the whole halo and 47 percent for the inner halo when 
we consider only haloes with no major mergers after formation. Our analysis suggests that 
spin flips (whose effects are not currently included in galaxy formation models) could be an 
important factor in the morphological transformation of disc galaxies. 
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1 INTRODUCTION 

The acquisition and evolution of angular momentum plays a cen- 
tral role in the formation and evolution of cosmic structure. Early 
work on the acquisition of angular momentum by virialized mat - 
ter clumps in a cosmological context dates back to iHovld (fl95lh. 
and t h e development of the li near tidal torque theory ( PeeblesI 
19691: jporoshkevichl Il97 0a"b': 'White' 'l984': 'Catelan & TheunsI 
1996al lbr see also IPorciani. Dekel & Hoffman 2002, and Schafeij 



yyell- r esolved objects 
20051 : IShawetalJ |200i 



;ts (e.g.J 
l2006l: Ih 



Bullock et al."200l': 'Avila-Reese et alj 
Hahnetal. 2007a b: BeU et al. 2001 



2010l: iMaccio et al.ll2007l: iMaccio. Dutton & van den BoscHUoOS 
Knebe & PoweJl2008l : lMunoz-Cuartas et alj201lh . 



These studies have usually focused on the distribution and 
evolution of the angular momentum magnitude. In contrast, the an- 
gular momentum vector direction is relatively less well studied, and 



20091) . This approach starts to bre ak down as structure growth 



becomes non- linear JWhitd Il984l) . and the subsequent evolu- 
tion of dark matter halo angular momentum is usually stud- 
ied using A'-body simulations. This subject too has a long 
and rich research history, with simulations improving in size 



(e.g. IWarren et al.l 


19921: iBailin & Steinmetz l2005l: lAlleood et alj 


20061 Shawetal 


20061: Havashi, Navarro & Sorineell 20071 


Bett et alj l2007l 


201of), or with other structures on differ- 


ent scales, such as galaxies (e.g. 'van den Bosch et al.' '200j; 



(e.g. PeeblesI 1971: 


Efstathiou & JonesI 19791: Davis et al 


1985 


Barnes & EfstathioiJ 


19871: Frenk et al. fl988l: IWarren et al. 


1992 


Cole &Lacevll 19961) 


Recent simulations have established the dis- 
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tribution of the angular momentum of dark matter haloes, and 
its evolution, extremely accurately, from very large numbers of 
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evolution of the mass in z = haloes, and showed that their 
spin direction changes due to non-linear evolution, with both the 
average deviation from the initial direction, and the scatter in that 
angle, increasing with time. 

The motivation for studying the angular momentum of haloes 
is the influence it is believed to have on the formation and evolu- 
tion of galaxies. In today's cosmological paradigm, in which the 
dark matter is a cold coUisionless particle (CDM), haloes grow 
hierarchically, through a series of mergers of ever larger objects. 
In the basic two-stage pictu re of galaxy formation dWhite & Reed 
1 19781 : [white & Frenkll99l|) , galaxies form and evolve within these 
haloes, with the pattern of simple hierarchical growth modified by 
the more complex physical processes available to the baryons as 
they cycle between gas and stars. It is usually assumed that the 
gas and dark matter have the same initial distribution, and thus the 
gas in a halo initially has the same angular momentum as the halo 
itself. The gas then collapses to form a (rotationally supported) 
disc galaxy, conserving its angular momentum. Thus, the size of 
the galactic disc is directly related to th e dark matter halo's angu- 
lar momentum jFall & Efstathiou|[T980l : IMo. Mao & Whitdl 19981 ; 
IZavala. Okamoto & Frenkll2008h . This basic picture is frequently 
impl emented in so-called 'semi-analytic' models of galaxy forma- 
tion ( IWhite&Frenkll99lh . in which the modelling of the baryonic 
processes is grafted onto the merger histories of dark matter haloes, 
derived either from an A'-body simulation or const ructed analyti- 
call y. This approac h has been recently reviewed by ISaughl ( l2006l) 
and lBensonloOlOh. and a comparison of different models has been 



and lcensonl ( 1/UtUI). and a compariso n ot airterent models nas been 
carried out by De Lucia et al.l ( |2010|) . It is important to emphasise 



that while these models incorporate the size of the dark matter an- 
gular momentum vector, they make no reference to its direction. 

It has long been know n that tidal forces can ind uce morpho- 
logical changes in galaxies jToomre & Toomr3ll972h . If the gravi- 
tational potential varies significantly over a short timescale, a galac- 
tic disc can be disrupted completely. Galaxy formation models 
thus assume that a sufficiently big galaxy merger event will de- 
stroy adisc,randomi^mgfi^^ 

(e.g.lToomrelll977l:lBameslll988l 1 1 9921 : iBames & Hemquistlll996l ; 
lHemauisllll99l Il993h . Indeed, this has been shown to occur 



in merger simulations of individual obje cts (e.g. |Na ab & Burkerj 
l2003l : lBoumaud. Jog & Co mbes 2005; Cox et alj2006 . 2008). The 
outcome of a merger depends on the gas richness of the pa r- 
ticipa nts (e.g. IStewart etal.ll2008l [200 % iHopkins et al. I l2009al lbl 
|20IO|) . and on the details of the star formation and feedback 
processes triggered by the merge r (e.g. Okamoto et alj l2005t 
[Zavala, Okamoto & Frenk 2008; Sc annapieco e t al. 2009). 

In this paper, we consider the evolution of the direction of the 
angular momentum vector (hereafter spin direction, for brevity) of 
dark matter haloes, a process that can affect the stability of a disc 
within the halo. Sudden, large changes in the halo spin direction 
are indicative of a significant disturbance to the halo. Such changes 
would usually accompany a halo mergei[|, which, in turn, could 
result in a galaxy merger within the halo, and potentially the de- 
struction of an existing galactic disc (depending on the details of 



' Such a spheroid is often distinguished from a so-called pseudobulge (e.g. 
iKoi-mendv & Kennicutj|2004l : lFreemaj2008l) . Pseudobulges are thought to 
form through secular evolution of the disc. 

^ In the case of a major merger, the result ing spin direction is corr elated 
with the net orbital spin of the progenitors jFaltenbacher et aOllOOSl) . Fur- 
themiore, there is a degree of coirelation in the infall directio ns of satel- 
lite haloes (e.g. falling along filaments. iKnebe et al. 2004. ; .Libeskind et alj 
l2005l ; lLovell et al.l201ll:lLibeskind et al.l201 ll) . 



the baryonic physics). However, it is also possible for tidal forces 
to disturb the halo without there being an immediate merger - for 
example, due to the flyby of a neighbouring halo. Recent work by 
ISinha & Hollev-B ockelmann (2 0110 has shown that halo flybys in- 
deed occur sufficiently frequently that they have a significant dy- 
namical effect on halo systems. In such a situation, the spin di- 
rection could change significantly, even if its magnitude does not. 
Such "spin flips" could have major consequences for the survival 
of a disc. The disturbance in the internal structure of the halo itself 
could lead to the distruction of the disc, as in a merger. However, it 
could also tor que the disc and chang e its spin direction without dis- 
rupting it (Os triker & Binnevlll989l) . causing it instead to become 
misaligned relative to the direction of new infalling material. In due 
course, the accretion of misaligned material could lead to the dis- 
ruption of the disc and the formation of a spheroid. Such spin flips 
provide a mechanism of spheroid formation that is not currently 
considered in galaxy formation models. 

A'-body and hydrodynamical simulations have shown that 
large, rapid changes to dark matter halo spin directions do indeed 
occur. IOkamotoetai] ( l2005h remarked that their simulated galaxy, 
which had formed a small disc, flipped its orientation. It then be- 
gan to accrete gas in a direction nearly perpendicular to the original 
disc which was sub sequently transformed into a bulge with a new 
disc later forming iRomano-Diaz et al] j2009l) analysed haloes in 
simulations both with and without baryons, at very high time res- 
olution. They found that, although the spin magnitude changes by 
a factor of ~ 2-3 after a major merger, the orientation can change 
much more drastically, by > 180°. Furthermore, such large changes 
in angular momentum orientation are not restricted to major merg- 
ers. In the simulations with baryons, the authors also found that the 
DM halo, stellar disc, and gas component can often flip orienta- 
tion with respect to each other as the sy stem evolves, even at late 
times where there are few major mergers. IScannapieco et al.](l2009l) 
found that misalignment of a stellar disc with the accreting cold gas 
can sometimes cause mass to transfer from the disc to a spheroidal 
component, sometimes destroying the disc (and sometimes allow- 
ing a new disc to be formed later). 

In this and subsequent papers, we use a dark matter A'-body 
simulation to assess the frequency of spin flip events occurring in 
the lifetime of haloes. We do not model baryonic physics, and in- 
stead concentrate on quantifying the amplitude and frequency of 
spin flips. Our aim in this first paper is to make an initial assess- 
ment of the importance of spin flips as a potential mechanism for 
the disruption of discs and the formation of spheroids. We focus 
here on those haloes whose mass at z = is similar to that of the 
Milky Way (MW). We perform a more in-de pth study on the dis- 
tribution of spin fl ips in a subsequent paper (iBett & Frenkll20I ll 
hereafter IPaper III) . 

The outline of this paper is as follows. In section |2] we de- 
scribe the A'-body simulation we use and our analysis procedure, 
including details of halo identification, merger trees, halo selec- 
tion, and the definition of the quantities of interest here: the frac- 
tional mass change, and the spin orientation change. We present 
our results in section |3] describing both the joint and cumulative 
distributions of spin flip and merger events, and investigate the fre- 
quency of spin flips over the course of halo lifetimes. We discuss 
our conclusions in section |4] 



' A prelimin ary analysis of this system, in the spirit of the present paper, 
can be seen in lBetil)2010l) . 
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^box ^ part 

A-'Mpc 
100 729 X 10^ 



lO'/r'M© r'kpc 



9.5ie 



2.4 



Table 1. Simulation parameters for the hMS simulation: box size, numbers 
and masses of particles, and gravitational softening 



0.75 



0.25 0.045 0.73 



1.0 



0"8 

0.9 



Table 2. Cosmological parameters (at ; = 0) for the hMS simulation used 
in this paper: cosmological density parameters H/o, the Hubble parameter 
h, the spectral index n, and erg is the linear theory mass variance in spheres 
of radius 8 /i" ' Mpc at z. = 0. As with the Millennium Simulation, these 
parameters wer e chosen to be good matches to the results of the 2dF galaxy 
redshift survey jColless et al...2001; Percival et al.„20(3)^ and the first yea r 
results of the WMAP microwave background satellite jSpergel et al.l20O3l) . 



2 SIMULATION DATA AND ANALYSIS 

In this section we describe the A'-body dark matter simulation we 
use and the associated halo catalogues and merger trees constructed 
to link each halo with its descendent in a subsequent output time. 
We calculate various halo properties, and use the merger trees to 
describe how these properties evolve over the lifetime of each halo. 



2.1 The hMS simulation, haloes and merger trees 

We use the hMS cosmological dark matter simulation. This was 
canied out using the same L-GADGET-2 code and ACDM cosmo- 
logica l parameters as the Millennium Simulation dSpringel et alj 
|2005|) . but with a smaller box size and higher resolutior0. Rel- 
evant simulation parameters are shown in Table [T] while the as- 
sumed cosmological parameters are shown in Table|2] Throughout, 
we refer to cosmological density parameters f2,(z) = p,(z)/Pc(z), 
in terms of the mass densit}|f] of component ; and the critical den- 
sity Pc(z) = 3H^{z)/(SnG), where H{z) is the Hubble parameter. 
We use a subscript zero to denote parameters evaluated at z = 0, 
and parameterise the present day value of the Hubble parameter as 
Ho = lOO/ikms-'Mpc"'. 

As with the Millennium Simulation, at each simulation snap- 
shot particle groups were identified on-the-fly according to the 
friends-of-frien ds algorithm (FoF ), with a linking length parame- 
ter of b = 0.2 l lDavis et"^l 19851) . Subsequently, self-bound sub- 
structu res within these grou ps were found using the Subfind algo- 
rithm dSpringel et all 1200 lb . Finally, the progenitors and descen- 
dents of each particle group were found, creating a 'merger tree' 
structure allowing haloes to be tracked over time. The merger tree 
algorithm (and asso ciated halo definition) used is that described in 
iHarker et all ilOOdt) . which was originally designed for use with 
the Galform semi-an alytic galaxy formation model and t he Mil- 
lennium SimulatiorQ ( iHellv et alj|2003l : lBower et al.ll2006h . Since 



The hMS simulation was previous l y used bvlNeto et al.l J2OO7I) .. 
2008l).lBovlan-Kolchin etal] l200Sl) . lLibeskind et alj<2009h . and 



201c ). 



Gao et al 



Bett et al. 



The equivalent mass density of the cosmological constant. A, is pA = 
Ac-/(S!tG). 

^ In particular, t hey correspond to the DHalo tables in the Millennium Sim- 
ulation database iLemson & the Virgo Consortiumll2006l) 



the construction of haloes and mergers trees from particle groups 
in discrete snapshots is essential for the current work, we will now 
describe this process in more detail. 

The preliminary set of haloes consists, at each snapshot, of 
the FoF particle groups, which in turn contain the self-bound sub- 
structures identified by Subfind (these include the main body of the 
halo, plus less massive subhaloes), as well as so-called "fuzz" parti- 
cles not gravitationally bound to any structure in the halo. It is well 
known that the purely spatial nature of the FoF algorithm allows 
multiple objects to be linked together spuriously, in the sense that 
although they are close together they might not necessarily be phys- 
ically connected. A common example is that of two close objects 
linked with a tenuous bridge of particles, which t he FoF algorithms 
identifies as a single "halo" (see lBett et al ] |2007l for a detailed dis- 
cussion and comparison of the effect of groupfinders on the mea- 
sured halo angular momentum and related properties). 

Including additional physical information - e.g. gravitational 
binding from Subfind, and temporal evolution from merger trees - 
allows the operational definition of a halo to be refined, to match 
better our physical intuition of what a halo is, and thus, when used 
in conjunction with semi-analytic models of galaxy formation, to 
allow better comparison s with both observation s and hydrodynamic 
simulations. Following IWechsler et al.l ( |2002|) . a "splitting" algo- 
rithm is applied to the basic FoF halo catalogues, whereby spuri- 
ously linked subhaloes are split off from their original FoF groups 
and identified as separate haloes in their own right. A subhalo is 
split off from its original FoF parent if it satisfies at least one of the 
following conditions: (1) The distance between the subhalo centre 
and the parent centre is more than twice the half-mass radius of the 
parent; or (2) the subhalo still has more than 75% of the mass it had 
when it was last identified as a separate halo. This yields halo cat- 
alogues containing more objects than the corresponding FoF cat- 
alogues. Ilettetal] 12003) showed that these 'merger tree haloes' 
are a great improvement on both the simple FoF groups and groups 
found from a simple spherical overdensity criterion. 

Merger trees for these haloes are constructed by tracking the 
particles that constitute the subhaloes between each snapshot, start- 
ing at early times and continuing to redshift z = 0. Each halo or 
subhalo can have at most one descendent in a later snapshot. The 
most bound 10% of a subhalo's mass (or 10 most bound particles if 
that is more massive) is located in the next snapshot. Occasionally, 
these particles might no longer reside in a subhalo in the next snap- 
shot: the subhalo might have temporarily dropped below Subfind's 
detection limit, or might be passing through a high density region 
and be interpreted as unbound matter around that density peak. In 
practice therefore, the next five snapshots are scanned to find the 
earliest time when these particles are again in subhaloes. The de- 
scendent subhalo is then identified as the subhalo containing the 
largest number of those most bound particles. The descendent of a 
halo is identified as the halo whose most massive substructure (i.e. 
the main self-bound halo component) is the descendent of its own 
most massive substructure. 

It is possible that a subhalo's mass ends up distributed between 
two (or more) subhaloes in a subsequent snapshot. While one will 
be identified as the descendent, the other will be left as a sepa- 
rate "orphan" object without a progenitor. This situation is known 
as a de-merger, and is a physical effect separate to the splitting of 
groups described above - here, sets of particles physically end up 
in separate objects as the simulation evolves. 

The end result of this process is a catalogue of haloes (groups 
of self-bound substructures) identified at each snapshot, with at 
most one descendent and one or more progenitors. Each halo iden- 
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tified at z = is the root of its own tree, whiclt branches into many 
progenitor haloes at preceding output times. In this paper, we study 
the evolution of properties of individual haloes that at z = have a 
mass corresponding roughly to that of the Milky Way. After iden- 
tifying an appropriate halo at z = 0, we track its evolution back by 
finding its most massive progenitor at the preceding snapshot, then 
finding the most massive of that halo's progenitors, and so on. 

It is important to note that, just like the halo definition and 
galaxy formation model, the halo merger tree algorithm is not 
by any means unique - even within a given A'-body simula- 
tion. The halo merger trees used for the " MPA" semi-analytic 
mo dels of the Millennium Si mulation (e.g. ISpringel et alj |2005| 
and IPe Lucia & Blaizoj 1 20071) in fact track the binding-energy- 
weighted mass in the Subfind subhaloes. Other methods that use 
splitting/stit ching algorithms similar to the one used here inc lude 
thosebve^Fakhouri & Ma ('2008','2009').'Ge nel et"ai] ( l2009h : see 
also iMaller et al. (2006). Tweed et al. (2009) provide a recent de- 
tailed study of halo definition and merger tree algorithms. 



2.2 Halo property catalogues 

Various properties of the haloes are computed at each output time. 
Properties are computed in the centre-of-momentum frame of each 
halo, and in physical rather than comoving coordinates. Each halo 
consists of a set of A'p particles, with each particle i having mass 
Mi = nip, position jc, and velocity v,. The halo mass is therefore 
Mh - m, - Nptrip. The halo centre is taken to be the location 
of the gravitational potential minimum of its most massive sub- 
structure, as found by Subfind. We define an approximate "virial" 
radius, ^vir, for the halc0 by growing a sphere from the halo centre 
and computing the density of the halo particles within. We locate 
/?vir at the radius at which the density enclosed drops below a cer- 
tain threshold value computed at that snapshot, ph = A;;(z)Pc(z). The 
threshold overdensity with respec t to critical, Ac(z), is dete rmined 
from the spherical co llapse model (pke^ Cole & Frenkll996h . using 
the fitting formula of ISrvan & Norman ( 1998b : 



A,(z) =lSn^ + 82 (Quiz) - 1) - 39 {Quiz) - 1)' 



(1) 



In the case of the flat ACDM universe used here, HmCz) = 
^Moo'^ /xiz) and pdz) = PcoXiz), where the expansion factor 
a = (1 -I- z)"' and we define ;^f(z) = t^ MpQ^ + ^ao for conveni ence. 

We compute halo energies as in lBett et alj ( I2007ll201ol) . The 
kinetic energy of a halo is given hy T - ^ Y^ili triiVj, and the poten- 
tial energy is computed as a double sum over a random sample of 
1000 particles, using the same smoothing kernel for gravitational 
softening as in the simulation itself. We only use the energies for 
some broad selection criteria (described below), and random sam- 
pling provides a good approximation. 

In this paper, we are mostly interested in the halo angular mo- 
mentum vector, / = iTiiXiXVi. We also define an inner halo an- 
gular momentum, /inner, using the particles within nnner = 0.25i?vii H 



2.3 Halo selection 

We need to select haloes at each snapshot from which reliable mea- 
surements of angular momentum can be made. The halo has to 



' Note that we do not use Ryii as a boundary for our halo. However, it 
provides a useful sc ale for the physica l halo size needed in other properties. 
^ This is similar to iBett et al.l fcOlOl) , although there we defined r-umQY — 
\0-° ^R„, 0.25i?vi, 



be well defined, that is, both well resolved (consisting of a suffi- 
ciently large number of particles), and reasonably relaxed (close 
to being virialised). Furthermore, the angular momentum magni- 
tude cannot be too small: if the angular momentum vectors of 
most particles are in opposite directions and cancel, then the net 
direction will be dominated by very few partic l es and will not 
be robust. We follow the approach of ISett et al. I ( l2010l) . defining 
a scaled angular momenturr|3 j as the ratio of the specific angu- 
lar momentum j to that of a single particle in a Keplerian orbit, 
/ = // \IGM\^R^-„, along with the analogous quantity for the inner 
halo Jinnei = dinner/ VGMinncr0.25/Jvir- As a basic Way of assesslng 
virialisation, we compute Q = 2T lU + 1, which should be around 
zero fo r a virialised halo . We use the same critical values for selec- 
tion as lSett et alj ( l201(il) : in particular, haloes that pass the follow- 
ing three criteria are retained: 



\Q\ 
logio J 



1000 

0.5 

-1.5 



(2) 
(3) 
(4) 



(When considering changes to the inner halo, the criteria for A'p 
and J are replaced by equivalent ones for A'p^inncr and Jinnei-) Note 
that these selection criteria are applied to haloes separately at each 
given snapshot, rather than once for their whole lifetime. Haloes 
can be excluded at one timestep (including at z = 0), but still re- 
tained for study at subsequent or preceding timesteps. 

In addition to the particle number cut above, we restrict our 
analysis to haloes that at z = have masses similar to that of the 
Milky Way halo. That is, we retain only haloes whose final mass is 
10506mp < Mo < 33 224rap, equivalent to the mass range 12.0 < 
logio(Mo/r'M0) < 12.5 (seeTable[B. 

A visual inspection of the coevolution of different halo proper- 
ties for individual haloes suggested that a further two selection cri- 
teria should also be applied. Firstly, the early life of a halo is very 
chaotic, with a high rate of mass accretion, mergers, and general 
instability in halo properties. So that our results are not dominated 
by this early period, before the halo has properly formed (in some 
sense), we restrict our analysis to the time period after the final time 
when M(z) < 0.5Mo (where Mo is the halo mass at z = 0). This cor- 
responds t o a commonly used s i mple definition of halo "formation" 
time (e.g.'Lacev & Cole" 1 993; 'Sheth & Tormen"2004'; 'Hahn et aU 
[2007a; Neto et al. 2007; Giocoli et al. 2007; Li, Mo & Gao 200a 
and references therein). Since stellar disks are unlikely to sur- 
vive t he early chaotic phase of halo formation dParrv. Eke & Fren3 
l2009l) . restricting attention to the period after the halo has formed is 
appropriate for investigating the frequency of spin flips that could 
alter the morphology of a disk galaxy. This condition must be borne 
in mind when interpreting the statistics that we present below. We 
plot the formation times of these haloes, according to this defini- 
tion, in Fig.[T] The peak in the distribution is around haloes forming 
atz a 1. 

Secondly, we found that there are some occasions in which 
the halo finder or merger tree algorithms make unphysical choices 
for which subhaloes to incorporate into which haloes. For exam- 
ple, if a satellite halo was orbiting near the edge of a halo, then 
it might "merge" at one snapshot, then be identified as a separate 
object again later, only to finally merge again afterwards. During 
such events, the angular momentum vector might appear to swing 
around wildly, as a large mass at large radius would be added to. 



^ Note that / is iden tical to the alternative spin parameter A' introduced by 
iBullock et all i200lh . modulo a factor of ^^2. 
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Redshift, z 
1 0.5 




cos 6{t) : 



J(tyJ{t-T) 

\J(t)\ \J(t-T)\ 



(6) 



Figure 1. Histogram of formation times (red) and initial detection ('start') 
times (blue) for haloes that have MW masses and pass our standard se- 
lection criteria at z = 0, in terms of lookback time, ?](,, and redshift. The 
formation and start times are computed using haloes at ; > that pass just 
the Np selection criterion. There is a histogram spike at each snapshot from 
z < 6.2; the blue spikes are offset slightly to make them visible. 



then removed from, the total halo /. Since such changes are not 
due to a physical change in the halo angular momentum, but in- 
stead are due to uncertainty in where to draw the halo boundary, 
we should exclude such events. It turns out that such events also 
cause large changes in the halo kinetic energy T, as the bulk ve- 
locity of the satellite halo will be incorporated into the main halo 
greatly increasing its net velocity dispersion. Thus, such events can 
be identified by considering the arithmetic change in the virialisa- 
tion parameter Q; since the potential energy U < and does not 
change much, an apparent sudden increase in T makes Q appear 
to decrease suddenly. By examining various cases, we chose to ex- 
clude eventj^that have Ag ^ -0.3. 

Finally, we note that we analyse the halo population over the 
redshift range z < 6.2; in any case, the effects we describe will 
be most visible at low redshift. As shown in Fig. [T] there are 46 
snapshots over this redshift range, with 23 over the period z < 1 . 



2.4 Evolution of halo properties 

Combining the merger tree data with the halo property catalogues 
at each snapshot, we can obtain the evolution of each halo property, 
for each halo identified at z = 0. We are most interested in the 
relationship between the change in halo mass and the change in 
halo spin orientation, from snapshot to snapshot. That is, we focus 
on two differential halo properties, the fractional mass change 

M(t) - M{t - t) 



Aiu(t) 



M(t) 



(5) 



and the angular change in spin orientation 



where t is the time at which the quantity is measured, and t is the 
timescale over which we measure the halo property change; the 
time t - T precedes the time t. 

In principle, we could simply look at the difference in halo 
properties at adjacent snapshot times. However, since the snapshots 
in the hMS are not evenly spaced in time, this would not be a fair 
way to analyse events in haloes at different times (the intersnapshot 
time varies between ~ 0.1-0.4 Gyr for z g 6). Instead, we choose a 
constant value for t, and simply linearly inteipolate the halo prop- 
erty in question between the values at the snapshots before and af- 
ter the time t - t. The simulation snapshots are in fact sufficiently 
closely spaced in time that this interpolation is accurate. 

We will refer to the property (or property change) of a given 
halo at a given snapshot as an event. We shall use some fiducial crit- 
ical values to divide the distribution of events to aid interpretation. 
We shall consider a spin direction change of at least 6*0 = 45° to be 
'large', and a fractional mass change of more than Ajig = 0.3 to cor- 
respond to a major mergei0 For the sake of brevity, we shall refer 
to events with AyU ^ 0.3 as minor mergers, even though they could 
be smooth accretion (i.e. not the merging with a satellite halo), or 
even mass loss. 

The choice of event timescale is non-trivial, since any charac- 
teristic halo dynamic timescale is likely to depend on halo mass and 
size (and therefore also on cosmology and time) - but we wish to 
use a single timescale for all haloes at all times, so that we can 
compare events at different times in different haloes fairly. The 
timescale we choose will therefore be only an approximation for 
the actual timescale of any particular halo. 

We consider the orbital timescale for a particle in a Keplerian 
orbit at the half-mass radius R1/2 of a model halo bounded by the 
radius enclosing the density Ac(z)pc(z). For a halo of a given mass, 
a concentration can be found by assu ming an NFW density profile 
l lNavarro. Frenk & White! 199611 1997 ) and using the redshif t depen- 
dent mass-concentration relation of Mufioz-Cuartas et^ ( I2OI l[) . 
The half-mass radius for a halo of a given concentration can then 
be found using the fitting formula of iLokas & MamonI ( 1200 ll) . al- 
lowing a timescale to be computed as 



Tl/2 



GM ■ 



(7) 



Since this effect is due to uncertainties in the halo boundary, we do not 
apply this exclusion criterion when considering the inner halo spin. 



For haloes in the mass range we consider in this paper, T1/2 varies 
from about 0.37 Gyr at z = 1 to about 0.63 Gyr at z = 0. (The 
values computed analytically using the fitting formulae outlined 
above agree with those measured directly from haloes in the sim- 
ulation.) We therefore take a fixed value for the event timescale of 
T = 0.5 Gyr, although we note that we do not expect our results to 
depend qualitatively on the exact value used; for some key results 
we show their depende nce on t. We will consider the choice of r in 
more detail in lPaper It is important to note that t is the timescale 
for our measurements of halo changes, and the physical timescale 
of flips or mergers can be much shorter. 



" A/j is only restricted to be < 1 ; a value of A/j = ^ means that the mass has 
increased by 50%. If A/i > j, then the halo has more than doubled in mass; 
we expect this to be rare, since we are, by definition, comparing with the 
most massive progenitor Negative values of A/i are possible, corresponding 
to mass loss between snapshots; A/j = -1 means that the halo has lost 50% 
of its previous mass. 
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Figure 2. Distribution of events as a function of fractional mass cliange, A^, 
and spin orientation change, cos d. Dotted lines mark the origin and dashed 
lines indicate our fiducial critical values for major mergers (A/.( 0.3) and 
large flips (B > 45°). 
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Figure 3. Cumulative distribution of events with spin misalignments of at 
least 00 degrees. The dashed line shows our fiducial value of 9o = 45°. 
We show results for a limiting merger fraction of A/io = 0.3 (red: major 
mergers; blue: minor mergers). The number of events in each case is written 
in the legend. 



Using the whole-halo selection criteria described in sec- 
tion 12.31 gives us a population of 35 279 events. When we select 
instead for the inner halo spin, we have 29 889 events. 



3 RESULTS 

3.1 The distribution of flips and mergers 

We start by examining how changes in the spin orientation of the 
halo correlate with changes in the halo mass that occur at the same 
time. 



3. 1. 1 Distribution of whole halo flip and merger events 

The distribution of events for MW-final-mass haloes is shown in 
Fig. (2] There are very few major mergers or large flips, with most 
events located around "no change" (cos 6 x: 1, A/j a 0). Most of the 
spread in cos 9 is located between no mass change and our fiducial 
threshold for major mergers (AyU = 0.3) 

Since we are interested in mergers and flips above and below 
some critical value, rather than at some value, it is useful to exam- 
ine the cumulative distribution functions (CDFs) of the data. The 
CDF of cos 6 is shown in Fig.[3] We can see that, if we consider just 
events without major mergers (AyU < 0.3), then only a very small 
fraction have large flips: about 0.7% have flips of 45° or more. If 
we select only major mergers, then since we have now excluded 
the main peak of the distribution we find a much higher proportion 
of events with large flips: about 17% have flips of at least 45° (al- 
though the major mergers themselves represent only 0.3% of the 
total event distribution). 

We can also consider the CDF of A/j (Fig.O. In this case, if 
we select just large flips, we find that the vast majority (93% of 
those with > 45°) coincide with minor mergers (A/j < 0.3). 
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Figure 4. Cumulative distribution of events with fractional mass change of 
A/Jo or less, for events with large spin flips (red: at least 45°; blue: at least 
90°). The dashed line shows our fiducial value of A/jq = 0.3, and the dotted 
line marks no mass change. 



3.1.2 Distribution of flips of the inner spin 

The distribution of flips of the inner halo angular momentum is 
more directly relevant when considering the stability of galaxies 
that might form within. The joint distribution of events as a function 
of the inner halo spin direction change and the mass change of the 
whole halo is shown in Fig. [5] In comparison to the distribution for 
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total halo spin flips, the inner halo exhibits a far greater spread to 
low-cos 9. 

Cumulative distributions are shown in the middle and right 
panels of Fig.[5] We find that the frequency of minor merger events 
(the blue line in the middle panel) that have a large inner spin flip 
is about 6.7%, which is a significant increase on that for total halo 
flips shown in Fig. [3] (0.7%). The fraction of major merger events 
that also have significant inner flips is slightly increased, to 26.6%. 
Selecting just large flips (right panel), the frequencies are similarly 
increased compared to the total halo flip distribution: 98.9% of flips 
of at least 45° coincide with minor mergers, dropping slightly to 
97.3% for flips of at least 90°. 

3.2 Spin flips over halo lifetimes 

While it is important to understand the overall frequency of flip 
events, and their tendency to correlate with mergers, we are also 
concerned with the frequency of spin orientation changes over the 
course of halo lifetimes. 



of haloes without major mergers since formation nevertheless have 
a flip of their total spin of at least 45° (185 haloes out of 2046). 
Since there are very few major mergers after fomation even for 
our total halo population (as shown in Fig.|2}, this figure does not 
change much if we do not apply the no-major-merger restriction: 
9.8% of such haloes have large flips, corresponding to 210 haloes 
out of 2146. On the other hand, 25% of the 100 haloes with major 
mergers since formation have spin flips of at least 45°. 

If we consider flips to the inner halo spin vector (right panel 
of Fig.[8l(, then as we have seen before, there is an increased likeli- 
hood for a halo to experience a significant flip. For haloes without 
major mergers after formation, 47% have large flips of their inner 
spin (946 out of 2006). Since large spin flips, particularly of the 
inner halo, are so common during the lifetimes of Milky Way-mass 
haloes, it seems unlikely that all such flips will result in the eventual 
destruction of a disc, although some fom of dynamical disturbance 
is to be expected. The elfect of the flips on the structure of the disc 
cannot, of course, be determined with our dark matter only simula- 
tions. 



3.2.1 Flips and coincident mergers as a function of flip duration 

An important question to answer is what is the likelihood of a halo 
exhibiting a spin flip (of a given magnitude 6o and measured over a 
timescale t) at some point during its lifetime? This can be further 
specified by restricting attention to spin flips that do (or do not) 
coincide with a major merger. 

We answer these questions in Fig. |6l for a range of values of 
6q and t, given our fiducial major merger threshold of AyUo = 0.3. 
As one would expect, the likelihood of getting a spin flip increases 
as one considers longer timescales or flips of smaller magnitudes. 
(The steps are an artefact of the discrete and irregular snapshot 
times, coupled with our interpolation scheme and relatively small 
halo population. Increasing t causes jumps as the snapshots used 
for interpolation change, and this occurs at different values of t for 
haloes in different snapshots. With a larger halo sample,, the lines 
become smooth curves, which we demonstrate in IPaper III ) 

Quantitatively, considering flips using our fiducial values of 
00 = 45° and r = 0.5 Gyr, we find that 10.5% of Milky- Way final 
mass haloes (172 haloes) experience such a flip at some point in 
their lives. If we consider just those for which such flips coincide 
with major mergers, this drops to just 0.9% (14 haloes). We find 
that 10.1% of haloes (166) experience a large flip without a major 
merger. 

We can construct a similar plot for changes to the inner halo 
angular momentum direction, which we show in Fig. [7] There is an 
increased tendency for haloes to have flips that do not coincide with 
major mergers, for all flip angles Oq and timescales t. For example, 
for our fiducial choice of t and 9o, we find that 58.5% (783 haloes) 
have large flips at some point in their lifetimes; a similar number 
have such flips without a major merger a (58.4%, 782 haloes). For 
those that have a major merger at the same time as such a flip, the 
fraction is still very low, at 1.1% (15 haloes). 

3.2.2 Flips in haloes without mergers 

Finally, we consider the particular case of haloes which have quiet 
merger histories, i.e. which do not have a major merger after their 
formation epoch (i.e. after z - 1 ; see Fig.lTJ. This case is interesting 
because it includes those haloes most likely to host a disc. 

The results are shown in Fig.|8](left panel). We find that 9.0% 



4 DISCUSSION AND CONCLUSIONS 

We have investigated the idea that spin flips - large and rapid 
changes in the orientation of the angular momentum of dark matter 
haloes - can occur without a major halo merger. These flips are a 
manifestation of strong tidal interactions and can be caused by mi- 
nor mergers or by flybys of a neighbouring object. Spin flips could, 
in principle, cause enough of a dynamical disturbance in the halo 
to disrupt or even destroy a galactic disc, perhaps resulting in the 
formation of a bulge or spheroid. Evidence for such dramatic out- 
comes have been seen in simulations o f the fomation of individ - 
ual galaxies (e.g. lOkamoto et al]|2005l ; IScannapieco et al]|2009l) . 
However, semi-analytic galaxy formation models do not take into 
account the potentially destructive effects of spin flips. The only 
processes that can transform discs into spheroids in current models 
are major mergers and disc instabilities caused by the accretion of 
matter onto the disc. 

Our goal in this paper has been to determine the frequency 
of spin flips during the lifetime of a galactic halo. We have dis- 
tinguished between flips affecting the entire halo and flips affecting 
only the inner parts of the halo which, at face value, would seem the 
most relevant for the stability of the disc. We have, for this initial 
exploration, chosen to focus on haloes that are roughly the mass of 
the Milky Way's halo at z = 0, i.e. between 10'^ and 10'^-' /i^'M© 
and that are reasonably re laxed, as would be expected for halos in 
which discs can form. In IPaper III we will extend this analysis to 
haloes of a much larger range in mass. 

We have found that, while the majority of what we have 
termed "events" (i.e. changes to a halo between a given snapshot 
tj and a preceding time /, - t) cause only small variations in both 
mass and spin direction, the distribution has a significant scatter 
and a large tail of significant variations in spin direction. The vast 
majority of large spin flips affecting the whole halo occur without 
an accompanying major merger (93 per cent of events with angu- 
lar change in spin direction 9 > 45° have a fractional mass change 
Ap < 0.3). However, such large halo-wide spin flips are a rare oc- 
cuiTence: only 0.7 per cent of non-major-merger events (Ap < 0.3) 
have f > 45°. 

Over the course of their lifetime (i.e. over the period after the 
halo has acquired half of its final mass), we find that 10.5 per cent 
of MW final-mass haloes experience at least one spin flip of 6 > 45° 



8 P. E. Bett (sr C. S. Frenk 



II, Number of mereer events per pixel 



180° 120° 



1.0 



0.5 



* 0.0 



-0.5 



-1.0 



■■ . "ft' 


1 


' ■■W^''. 

- ' 

, , , , 1 . , , , ^^^ ' . 





1.0 -0.5 



0.0 

An 



0.5 1.0 




Figure 5. Left: event distribution as a function of tlie inner halo spin flip versus tlie total lialo fractional mass change. Middle: Cumulative distribution function 
for merger events with inner spin misalignments of at least 6i„nei o degrees, for major merger events (red) and minor mergers (blue). Right: Cumulative 
distribution of events with mass changes < A/Jo, for spin flips of at least 45° (blue) and at least 90° (red). 
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Figure 6. The fraction of haloes with Milky-Way masses at z = that have at least one flip of at least 0o and duration r (left). Middle: same, but adding the 
constraint that the flip must coincide with a major merger > A/io, with A;jo = 0.3. Right: same, but the flip must not coincide with a major merger 



0.0 0.5 



2.0 00 0.5 rO 1.5 2.0 




0.0 0.5 



T/Gyr l/Gyr l/Gyr 

Figure 7. As Fig.|6] but for flips in the inner halo angular momentum. 
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Figure 8. The cumulative fraction of MW final-mass lialoes that have a spin flip (left: whole halo, right: inner halo) of at least 60 degrees at some point in 
their lifetime (black line). The other two lines correspond to the cases when the halo does (red) or does not (blue) have a major merger {Afi > 0.3) during its 
lifetime. (The black fine is almost coincident with the blue line.) 



with a timescale of t = 0.5 Gyr; 10.1 per cent of the haloes expe- 
rience such a flip without it coinciding with a major merger. These 
percentages increase for longer timescales and smaller minimum 
angle change. Finally, we find that 9 per cent of the haloes that 
have not had any major merger after formation nevertheless have at 
least one spin flip of 45° or more. 

The spin of the inner halo is subject to larger and more fre- 
quent changes in direction than the total halo spin, but like the total 
halo spin, inner spin flips also occur mainly without an accompa- 
nying major merger. Over half of the haloes have large inner spin 
flips at some point in their lifetimes without these coinciding with 
a major merger. For the haloes that do not experience any major 
mergers after formation, 47% experience a large inner halo spin 
flip. Large spin flips occur sufficiently frequently that they could 
have a significant impact on the evolution of the galactic baryonic 
material. 
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Our results suggest that a more complete understanding of the 
stability and resilience of galactic discs will require looking be- 
yond mergers (major or minor) and internal instabilites and should 
include the role of spin flips which, as we have seen, can be quite 
common for the inner halo. The survivability of discs is not just de- 
term ined by the halo potential. As has been sho wn in both models 
(e.g. 'Stewart et al."2009'; 'Hopkins et al."2009a"b^ and simulations 
(e.g. Okainoto et al. 2005; Scannapieco et al. 2009), the details of 
the baryonic physics - gas fraction, strength of supernova feedback, 
and other types of interaction between stars and gas - play a ma- 
jor role in whether a galactic disc lives or dies, or even can reform 
afterwards. Nevertheless, the behaviour of the underlying dark mat- 
ter plays a critical role in galaxy formation. To understand how spin 
flips influence the evolution of discs will require full baryon physics 
simulations at high resolution. The handful of examples of simula- 
tions we have mentioned here already demonstrate that this process 
is both important and tractable. 
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